% for a 64-channel EEG:
EEGnames = {
    'Fp1'    'Fpz'    'Fp2'    'F7'    'F3'    'Fz'    'F4'    'F8'    'FC5'    'FC1'    'FC2'    'FC6'    'M1'    'T7'    'C3'    'Cz'    'C4'    'T8'    'M2'    'CP5'    'CP1'    'CP2'...
    'CP6'    'P7'    'P3'    'Pz'    'P4'    'P8'    'POz'    'O1'    'O2'    'EOG'    'AF7'    'AF3'    'AF4'    'AF8'    'F5'    'F1'    'F2'    'F6'    'FC3'    'FCz'    'FC4'    'C5'...
    'C1'    'C2'    'C6'    'CP3'    'CP4'    'P5'    'P1'    'P2'    'P6'    'PO5'    'PO3'    'PO4'    'PO6'    'FT7'    'FT8'    'TP7'    'TP8'    'PO7'    'PO8'    'Oz'
    };

% The matrix is
Emap = [
    35 92  1 % Fp1
    50 95 2 % Fpz
    65 92  3 % Fp2
    17 79  4 % F7
    34 73  5 % F3
    50 69  6 % Fz
    66 73  7 % F4
    83 79  8 % F8
    20 62  9 % FC5
    40 59  10 % FC1
    60 59  11 % FC2
    80 62  12 % FC6
    0 34   13 % M1
    10 50  14 % T7
    29 50  15 % C3
    50 50  16 % Cz
    71 50  17 % C4
    90 50  18 % T8
    100 34 19 % M2
    20 37  20 % CP5
    39 40  21 % CP1
    61 40  22 % CP2
    80 37  23 % CP6
    17 22  24 % P7
    33 27  25 % P3
    50 30  26 % Pz
    67 27  27 % P4
    83 22  28 % P8
    50 15  29 % POz
    38 3   30 % O1
    62 3   31 % O2
    % EOG
    28 90  0 % AF7
    39 87  0 % AF3
    61 87  0 % AF4
    72 90  0 % AF8
    25 75  0 % F5
    42 70  0 % F1
    58 70  0 % F2
    75 75  0 % F6
    30 60  0 % FC3
    50 59  0 % FCz
    70 60  0 % FC4
    18 50  0 % C5
    39 50  0 % C1
    61 50  0 % C2
    82 50  0 % C6
    30 38  0 % CP3
    70 38  0 % CP4
    25 25  0 % P5
    42 28  0 % P1
    58 28  0 % P2
    75 25  0 % P6
    29 8   0 % PO5
    38 13  0 % PO3
    62 13  0 % PO4
    71 8   0 % PO6
    11 64  0 % FT7
    89 64  0 % FT8
    12 36  0 % TP7
    88 36  0 % TP8
    24 15  0 % PO7
    76 15  0 % PO8
    0 50   0 % Oz
    
    ];

% 

SPF = {
    [2 4 5 6]; %'1';
    [1 3 6]; %'2';
    [2 6 7 8]; %
    [1 5 9]; %
    [1 4 6 9 10]; %
    [5 7 10 11]; %
    [3 6 8 11 12]; %
    [3 7 12]; %
    [4 5 10 14 15]; %
    [5 6 9 11 15 16]; %
    [6 7 10 12 16 17]; %
    [7 8 11 17 18]; %
    [14 20 24]; %
    [13 9 15 20]; %
    [9 10 14 16 20 21]; %
    [10 11 15 17 21 22]; %
    [11 12 16 18 22 23]; %
    [12 17 19 23]; %
    [18 23 28]; %
    [14 15 21 24 25]; %
    [15 16 20 22 25 26]; %
    [16 17 21 23 26 27]; %
    [17 18 22 27 28]; %
    [20 25]; %
    [20 21 24 26 29]; %
    [21 22 25 27 29]; %
    [22 23 26 28 29]; %
    [23 27]; %
    [25 26 27 30 31]; %
    [29 31]; %
    [29 30]; %
    };

pairs = {
    [1 3];
    [2 2];
    [4 8];
    [5 7];
    [6 6];
    [9 12];
    [10 11];
    [13 19];
    [14 18];
    [15 17];
    [16 16];
    [20 23];
    [21 22];
    [24 28];
    [25 27];
    [26 26];
    [29 29];
    [30 31];};
    

% =====================
corrections = [1 2 3 4 8 14 18 24 28 30 31];
r = 43;
center_x = 50;
center_y = 48;
% for i=1:1:360
%     x(i) = r.*cos(2.*pi.*i./360)+center_x;
%     y(i) = r.*sin(2.*pi.*i./360)+center_y;
% end
% plot(x,y)

for i=corrections
    r_old = sqrt((Emap(i,1)-center_x).^2+(Emap(i,2)-center_y).^2);
    Emap(i,1) = (Emap(i,1)-center_x)./r_old.*r+center_x;
    Emap(i,2) = (Emap(i,2)-center_y)./r_old.*r+center_y;
end

